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ABSTRACT 

A comprehensive analysis of extended emission line region (EELR) around quasars is presented. New 
Subaru/Suprime-Cam observation is combined with literature search, resulting in a compilation of 81 
EELR measurements for type-1 and type-2 quasars with associated active galactic nucleus (AGN) and 
host galaxy properties. It is found that EELR phenomenon shows clear correlation with Eddington 
ratio, which links EELR to the constituents of the principal component 1 (PC 1), or eigenvector 1, 
of the AGN emission correlations. We also find that EELR is preferentially associated with gas-rich, 
massive blue galaxies. It supports the idea that the primary determinant of EELR creation is the gas 
availability and that the gas may be brought in by galaxy merger triggering the current star formation 
as well as AGN activity, and also gives an explanation for the fact that most luminous EELR is found 
around radio-loud sources with low Eddington ratio. By combining all the observations, it is suggested 
that EELR quasars occupy the massive blue corner of the green valley, the AGN realm, on the galaxy 
color - stellar mass diagram. Once a galaxy is pushed to this corner, activated AGN would create 
EELR by the energy injection into the interstellar gas and eventually blow it away, leading to star- 
formation quenching. The results presented here provide a piece of evidence for the presence of such 
AGN feedback process, which may be playing a leading role in the co-evolution of galaxies and central 
super-massive black holes. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: ISM — galaxies: star formation 
— quasars: emission lines 



1. INTRODUCTION 

Some active galactic nuclei (AGNs) are surrounded 
by massive ionized nebulae extending to several tens of 
kiloparsec from the central engines. The first system- 
atic study of such extended emission-line region (EELR) 
around quasars in combin ation with nuclear spe ctral 
properties was presented by IBoroson fc Okd (J1984D and 
iBoroson et al.l (|1985| ). They found that nearly half of 
their 24 objects are EE LR quasars (quasars with de- 
tectable EELR, following IFu fc Stocktonll2007D and that 
the EELR quasars have larger [O III] A5007 equivalent 
width (EW), broader and bumpier Balmer lines, and 
weaker Fe II emission compared to non-EELR quasars. 
In addition, the EELR quasars are predominantly as- 
sociated with steep-spectrum radio emission with ex- 
tended, lobe-dominant morphology, in contrast to the 
non-EELR quasars associated with flat-spectrum radio 
emission with compact, core-dominant morphology or 
without strong radio emission. 

The correlation between EELR phenomenon 
and radio emission was further confirmed by 
iStockton fc MacKentvl ()1987[ ). who obtained narrow- 
band [O III] line images of 47 quasars. Above the limiting 
luminosity of extended [O III] emission Le[oiii] = 5x 10^^ 
erg s~^, 10 out of 26 (10/26) of their steep-spectrum 
radio quasars (SSRQs) have EELR, while majority of the 
flat-spectrum radio quasars (FSRQs) and the radio-quiet 
quasars (RQQs; 7/7 of the FS RQs and 13/14 of th e 
RQQs) are non-EELR sources ()Fu fc StocktonI I2009D . 
They also found a correlation between nuclear [O III] 
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luminosity and the presence of EELR. iFu fc StocktonI 
(|2007f ) analyz ed the rest-frame ultraviolet s pectra of 12 
SSRQs in the IStockton fc MacKentvl ([19871 ) sample and 
found that N V A1240 emission is suppressed in EELR 
quasars, which they claimed points to lower metallicity 
of the broad-line region (BLR) compared to non-EELR 
quasars. They did not find any systematic difference in 
other AGN properties such as continuum luminosity and 
mass of super-massive black holes (SMBHs) between the 
EELR and non-EELR SSR Qs. 

On the other hand, Huse mann et al.l (|2008l ) observed 
20 quasars dominated by RQQs and found EELR around 
8 of them. The correlation between nuclear Fe II EW, 
H/3 line width, and EELR was confirmed in their sample 
too. In addition, they found a dependence of EELR on 
SMBH mass (A/bh) and reported that A/bh = 10^-^Mq 
separates EELR quasars with larger A/bh from non- 
EELR quasars with smaller A/bh- Se arches for EELR 
around R QQs were also presented by iHumphrev et al.l 
(,2010,) and iVillar-Martin et al.l (|2010l 120111) focusing on 
obscured, type-2 quasars. They found EELR around 10 
out of 20 objects and claimed that luminous EELR is not 
necessarily associated with powerful radio sources. 

Observational constraints on the emergence of EELR 
can possibly be a key to understand the co-evolution of 
galaxies and SMBHs. It is now widely accepted that 
the tight correlation between S MBH mass and bulge 
mass of their host g alaxies (e.g.. IMagorrian et al.lll998l : 
IHaring fc Rixl I2004D indicates the intimately-connected 
evolution of the two systems. The most compelling 
process of interaction between galaxies and SMBHs is 
the AGN feedback, in which the energy emitted by 
AGN regulates the mass accretion to SMBHs and drives 
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galactic wind that can expel cold gas from and quench 
star formation in host galaxies. The AGN feedback 
could also give a solution to the two major problems 
in the current galaxy formation models based on the 
A cold dark matter (CDM) theory: first is the "over- 
cooling" problem in which much more massive galaxies 
are formed in the numerical simulations than observed 
due to too efhcient gas cooling, and second is the "in- 
verted color - mass - morphology relation" problem in 
which massive galaxies arc predicted to be blue and 
disk-dominant ra ther than red and sph eroid-dominant 
as observed (e.g.. iSomerville et al.ll2008f ). However, ob- 
servational evidence is still lacking on whether the AGN 
feedback is actually present or not and how it works if 
present. While measurements of galaxy evolution im- 
printed on the s tellar-mass ( l uminosity) function (e.g 



Bell et al.l [200l iFaber et all [20071 iBrown et al.l 



2008f) especially at its mass i ve end (iConselice et alJ 
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its integration (the ext ragalactic background light; e.g, 
IMatsuoka et al.l l2011b[ ) , or analyses of the nature of 
the AGN p henomenon (e.g., Laor 2000; Boroson 2002; 
IMatsuoka e t al. 2007, 2008; Vestcrgaard et al. 2008) can 
put indirect constraints on the co-evolution, EELR pos- 
sibly provides more direct clues as to the interrela- 
tion between galaxies and SMBHs. The previous mea- 
surements indicate that EELR gas is largely ionized 
by AGN rather than stellar populatioii (iFu fc Stockton! 
[200l iHusemann eTalllmol : iKeel et al.l[2012t ). hence ^ 
may be witnessing the AGN feedback in action in this 
:alaxy -wide ionization process. Indeed, iGreene et al.l 
20 lit) recently showed that AGN photoionizes and kine- 
matically disturbs the interstellar medium throughout 
the entire host galaxies of luminous type-2 quasars. 

In this paper, new observation of EELR around five 
quasars at z '-~^ 0.3 is first presented. Then it is combined 
with the previous observations of EELR and associated 
AGN and host galaxy properties obta ined mainly frorn 
the Sloan Digital Sky Survey (SDSS; I York et al.ll2000l) 
archive. By using this compiled sample as well as the 
control samples of SDSS quasars and galaxies, we aim to 
perform a comprehensive analysis of EELR in the context 
of galaxy evolution. A cosmology with Hq = 70 km s~^ 
Mpc"-^, flu = 0.3, and ^a = 0.7 is assumed throughout 
this work. Luminosity and mass are given in units of erg 
s~^ and the solar mass (M©) unless otherwise noted. All 
the presented magnitudes are on the AB system. 

2. DATA 

2.1. New Observation 

The observed quasars were selected from the SDSS 
Data Release 7 (pR7) quasar catalog (jSchneider et al.l 
120101 : iShen et all 120111 1. The simple selection criteria 
were set based on redshift (z ~ 0.3 in order that [O III] 
A5007 falls in the narrow-band filter transmission de- 
scribed below) and nuclear [O III] luminosity (iN[oiii] > 
j^q42 gj.g s-i)_ Among the objects meeting the criteria, 
those with various SMBH mass and radio loudness and 
with the highest [O III] luminosity were picked up. The 
properties of the observed quasars are summarized in Ta- 
ble [TJ They effectively complement the existing measure- 
ments on the AGN principal-component plane explored 
below. 



The targets were obser ved with the prime-fo cus imag- 
ing facility Suprime-Cam (jMivazaki et al.ll2"002D mounted 
on Subaru telescopo The instrument covers a wide 
field-of-view of 34 x 27 arcmin^ with ten 2k x 4k CCDs, 
with a pixel scale of 0.20 arcsec. The narrow-band NAgse 
(effective wavelength Aeff = 6570 A, effective bandpass 
AAoff = 140 A) and the broad-band Re (XcS = 6530 
A, AAoff = 1110 A) filters were used to trace the red- 
shifted [O III] line and underlying continuum, respec- 
tively. The observation was carried out on 2011 May 3-4 
under the open-use program SllA-028. The sky condi- 
tion was sometimes non-photometric, but the measured 
flux of celestial objects was mostly stable within 10 % 
variation. The seeing was 0.8-1.1 arcsec. For each tar- 
get, nine 300-sec or three 900-sec exposures in NAgse 
and nine 100-sec or three 300-sec exposures in Re were 
obtained, depending on the source brightness. In addi- 
tion, a spectrophotometric standard star SA 104-335 was 
observed with the same instrument configuration as the 
target observation. 

Data reduction was performed with the Suprime -cam 
Deep Field Reduction (SDFRED: [Ouchi et al.ll200l soft- 
ware in a standard manner, including bias subtraction, 
flat flelding, distortion correction, sky subtraction, mask- 
ing bad pixels, and stacking. Photometric calibration 
in NAgse was achieved by referring to the observed flux 
of SA 104-335, whose intrinsic magnitude was derived 
by convolving its spectrophotometry data (jStond [T996h 
with the fllter transmission function. Then the Re im- 
ages were calibrated relative to NAgse by requiring the 
mean (i?c — NAgse) color to be zero for galaxies detected 
with high signal-to-noise ratios (point sources were not 
used because stellar Ha absorption could affect this cal- 
ibration). The scatter of ^0.05 mag was found around 
the mean between Re and NAgse magnitudes in the cal- 
ibration. 

The reduction process then proceeded to the contin- 
uum subtraction from the NAgse images. For this pur- 
pose, the nuclear spectra of the quasars were retrieved 
from the SDSS DR7 archive! (an example is shown in 
Figure [1]) . They were measured within SDSS spectro- 
scopic fibers of 3-arcsec aperture. For each spectrum, 

^ This work is based in part on data collected at Subaru Tele- 
scope, which is operated by the National Astronomical Observatory 
of Japan (NAOJ). 
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versity, the Joint Institute for Nuclear Astrophysics, the Kavli In- 
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entist Group, the Chinese Academy of Sciences (LAMOST), Los 
Alamos National Laboratory, the Max-Planck-Institute for Astron- 
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EELR and the AGN Feedback 



TABLE 1 

Targets of Subaru/Suprime-Cam observation. 



r-band 
mag.'' 



Rodshift 



log Lf^lO III] 
(erg s-l) 



log Mbh 

(Me) 



Rb 



log Le[o III] 
(erg s-l) 



SD9S 091401.75+050750.6 


17.32 


0.301 


42.85 ± 0.01 


9.36 ± 0.10 


L 


42.17 ± 0.08 


SDSS 095456.89+092955.7 


17.92 


0.298 


42.93 ± 0.01 


8.56 ± 0.08 


L 


< 41.47 


SDSS 113949.47+402048.5 


17.52 


0.314 


42.97 ± 0.01 


8.21 ± 0.08 


Q 


< 41.48 


SDSS 150740.92+445331.5 


18.16 


0.314 


42.69 ± 0.01 


7.34 ± 0.14 


L 


< 41.53 


SDSS 150752.66+133844.5 


17.79 


0.322 


42.53 ± 0.01 


9.67 ± 0.02 


Q 


41.83 ± 0.08 
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Note — ^SDSS PSF magnitude. The measurement error is < 0.02 mag. ''Radio loudness. Q; radio quiet, L: radio loud. 

the continuum levels beneath [O III] A4959 and A5007 
were estimated by interpolation between continuum win- 
dows at both (shorter- and longer-wavelength) sides of 
the lines, and then the spectrum was decomposed into 
[O III] and continuum. Convolving them with the fil- 
ter transmissions gives us the relative fluxes of [O III] 

and continuum in NAgse (/nAbsI ^'^^ ^^K^^^ and in R^ 
(/jj ' and /^°"*) for each object. Using these measures, 
accurate continuum subtraction was achieved as follows. 
The radial profiles of all the quasar nuclei on the obtained 
images are consistent with the point spread functions 
(PSFs) derived from nearby stars, hence they are domi- 
nated by the radiation from central unresolved sources. 
The NA656 and i?c images were decomposed into nuclear 
(NAg^^, i?™'=) and extended (NA^^^, Rf^^) components 
by the best-fit Moffat profiles, fitted within a 10-pixel 
(2-arcsec) aperture on top of the quasar peak positions. 
Then the nuclear [O III] distributions are given by 



NA-^ X 



„[OIII] 
•'NAese 



j^[OIII] . /.cont ' 
/NAese "*" -/NAese 



while the extended [O III] distributions are 
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The second term in the expression ([5|), which is equiv- 
alent in principle to the ratio AAeff,R„ / (AAeff,R^ — 
AAcff,NA656)j corrects for the continuum oversubtraction 
of the first term due to the [O III] themselves contained 
in i?™*. Note that the expression ([1]), although it is most 
accurate, is not applicable to the extended components 
since the relative fluxes of [O III] and continuum are 

different from /jtj^ and /^X* measured in the SDSS 
spectroscopic fibers. It is worth re-emphasizing that the 
spectral information is quite useful in the above process 
as an accurate cross-calibrator of the images taken in the 
different filters. 

2.2. Data Compilation 

Thanks to the massive spectroscopic observations 
conducted by the SDSS, nuclear properties are now 
known for many quasars with EELR measurements 
in the literature. We compile the EELR measure - 
ments of typ e-1 quas a rs fro m iBorospn et al.1 (Il985f). 

tockton & M acKentvl (|1987f ). and iHusemann et all 
008>,2010.). The continuum an d line luminosities given 
in lStockton fc MacKentvl (|1987f ) are converted to the val- 
ues in the cosmology adopted here. The measurements of 
[Husemann ct al., (,2008, ) are corrected in accordance with 
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Fig. 1.— Nuclear spectrum of SDSS 1139+4020 (solid line) with 
the arbitrary-scaled transmission functions of broad-band Re and 
narrow- band NAese filters (dotted lines). 

their re-analysis (Husemann et al., in preparation). Then 
we cros s-match the quas ars with the SDSS DR7 quasar 
catalog (|Shen et al.ll201lD and obtain the nuclear proper- 
ties such as continuum and line luminosity, SMBH mass, 
and radio characteristics based on the Faint Images of the 
Radi o Sky at Twenty cm survey (FIRST; IBecker et al.l 
119951 ). For some objects without counterparts in the 
SDSS catalog, the nuclea r properties are t aken from the 
above papers as well as iFu fc: StocktonI ()2007f ). Com- 
bined with the new observation presented in this paper, 
we end up with 61 type-1 quasars with EELR measure- 
ments (27 detections and 34 non-detections) and with 
a variety of nuclear properties. The compiled sample is 
listed in Tabled We assign radio flag to RQQs, flag 1 
to fiat-spectrum (radio spectral index a^ < 0.5) or core- 
dom inant (based on the FIRST morphologv: lJiang et al.l 
l2007f l radio quasars, fiag 2 to steep-spectrum {an > 0.5) 
or lobe-dominant radio quasars. The SDSS quasars are 
separated into radio-quiet and ra dio-loud sources a t the 
radio- loudness parameter R = 5 (jShen et al.ll201ll ). For 
simplicity the objects with radio flag 0, 1, 2 are hereafter 
referred to as RQQs, FSRQs, and SSRQs, respectively, 
since there is a well-known, strong correlation between 
the radio morphology and spectral index. 

We also collect the EELR measurements of 20 type-2 
qua sars (10 detection s and 1 non-d etections) presented 
by 'Hump hrev et all (|2010[ ) and iVillar-Martin et al.l 
([2010, 2011]). Since they w ere selected froin the S DSS 
type-2 quasars identified by iZakamska et al.l ()2003| ) , lu- 
minosity and EW of emission lines arising from the 
narrow-line region (NLR) as well as optical ugriz mag- 
nitudes of the host galaxies are available in the archive. 
We derive the rest-frame g — i color and stellar mass of 
the host galaxies from the SDSS model magnitudes and 
spectros copic redshifts, using a fc-correction code devel- 
oped bv lBlanton fc Row cis (2007). The contributions of 
relatively strong [O II] A3727 and [O III] A4959, A5007 
lines from quasars are removed before the calculation by 
subtracting Amx = -2.5 log (1 -I- fx,Y [EWy/wjc]) 
from X (u, g, r, i, or z)-band magnitude if the line Y is 
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TABLE 2 
Nuclear properties of type-1 quasars with EELR measurements. 







log ALasioo 


log Mbh 




log ■f'N[OIII] 


log -^E[OIII] 




Quasar 


Redshift 


(ergs-l) 


{Mq) 


R^' 


(erg s-i) 


(erg s-l) 


Ref. 


4C 15.01 


0.450 


45.61 


9.24 


2 


43.58 


< 42.56 


2 


PG 0026+129 


0.142 


44.96 ± 0.11 


8.98 





42.53 ± 0.09 


41.80 ± 0.14 


3 


PG 0050+124 


0.061 


44.31 ± 0.11 


8.00 





42.26 ± 0.09 


< 43.0 


3 


PG 0052+251 


0.155 


44.98 ± 0.10 


8.55 





42.70 ± 0.09 


41.75 ± 0.13 


1, 2, 3 


SDSS 0057+1446 


0.172 


44.93 ± 0.01 


9.38 ± 0.02 





42.23 ± 0.03 


< 43.0 


3 


HE 0132-0441 


0.154 


44.69 ± 0.10 


7.73 





42.15 ± 0.16 


< 43.0 


3 


NAB 0137-01 


0.334 


45.17 ± 0.01 


9.01 ± 0.02 





42.62 ± 0.02 


< 41.53 


2 


SDSS 0155-0857 


0.165 


44.43 ± 0.02 


8.82 ± 0.04 





42.15 ± 0.02 


41.43 ± 0.19 


3 


HE 0157-0406 


0.218 


44.72 ± 0.09 


8.52 





42.09 ± 0.10 


< 43.0 


3 


Mrk 1014 


0.163 


44.78 ± 0.01 


8.06 ± 0.05 


1 


42.74 ± 0.02 


42.22 ± 0.14 


2, 3 


OI 287 


0.444 


45.08 ± 0.01 


9.09 ± 0.21 


2 


43.51 ± 0.01 


< 42.47 


2 


SDSS 0836+4426 


0.254 


45.29 ± 0.01 


8.73 ± 0.06 





43.28 ± 0.01 


43.08 ± 0.09 


3 


SDSS 0914+0507 


0.301 


44.84 ± 0.01 


9.36 ± 0.10 


2 


42.85 ± 0.01 


42.30 ± 0.08 


5 


SDSS 0948+4335 


0.226 


44.76 ± 0.01 


8.43 ± 0.03 





42.44 ± 0.01 


< 43.0 


3 


SDSS 0954+0929 


0.298 


44.54 ± 0.01 


8.56 ± 0.08 


2 


42.93 ± 0.01 


< 41.60 


5 


Ton 28 


0.329 


45.44 ± 0.01 


8.03 ± 0.43 





42.42 ± 0.02 


< 42.34 


2 


4C 13.41 


0.241 


45.41 ± 0.01 


9.54 ± 0.05 


2 


42.40 ± 0.02 


< 41.41 


2 


HE 1029-1401 


0.086 


44.97 ± 0.05 


8.70 







41.78 ± 0.12 


4 


PG 1049-005 


0.359 


45.59 ± 0.01 


9.17 ± 0.03 





43.47 ± 0.01 


< 42.31 


2 


4C 61.20 


0.421 


45.26 ± 0.01 


9.10 ± 0.09 


1 


43.34 ± 0.01 


< 42.55 


2 


4C 10.30 


0.422 


44.78 ± 0.01 


9.33 ± 0.05 


2 


42.75 ± 0.01 


< 41.62 


2 


3C 249.1 


0.313 


45.42 


8.96 


2 


43.46 


42.85 ± 0.02 


1, 2 


SDSS 1131+2632 


0.244 


44.56 ± 0.02 


9.08 ± 0.04 





42.81 ± 0.01 


42.11 ± 0.14 


3 


SDSS 1139+4020 


0.314 


44.63 ± 0.01 


8.21 ± 0.08 





42.97 ± 0.01 


< 41.61 


5 


4C 49.22 


0.334 


44.70 ± 0.01 


8.50 ± 0.01 


1 


42.45 ± 0.01 


< 41.81 


2 


GQ Comae 


0.165 


44.15 ± 0.01 


8.36 ± 0.03 





42.58 ± 0.01 


40.19 ± 0.30 


1, 2 


PKS 1217+023 


0.240 


45.13 ± 0.01 


8.87 ± 0.08 


2 


42.78 ± 0.01 


41.82 ± 0.16 


2, 3 


4C 25.40 


0.268 


44.73 ± 0.01 


8.86 ± 0.03 


2 


42.66 ± 0.01 


42.61 ± 0.02 


2 


30 273 


0.158 


46.04 ± 0.09 


9.06 


1 


42.87 ± 0.10 


< 41.73 


1, 2, 3 


SDSS 1230+6621 


0.184 


44.53 ± 0.01 


8.26 ± 0.03 





42.66 ± 0.01 


41.90 ± 0.15 


3 


HE 1228+0131 


0.117 


44.95 ± 0.09 


8.29 





42.23 ± 0.11 


< 43.0 


3 


SDSS 1230+1100 


0.236 


44.73 ± 0.01 


8.60 ± 0.05 





42.64 ± 0.01 


< 43.0 


3 


PG 1307+085 


0.154 


44.80 ± 0.01 


8.59 ± 0.05 





42.69 ± 0.01 


< 41.06 


2 


B2 1425+267 


0.364 


45.14 ± 0.02 


9.66 ± 0.03 


2 


42.99 ± 0.01 


42.37 ± 0.04 


1, 2 


PG 1427+480 


0.221 


44.62 ± 0.01 


7.95 ± 0.03 





42.52 ± 0.01 


< 43.0 


3 


SDSS 1444+0633 


0.208 


44.48 ± 0.01 


8.33 ± 0.03 





42.44 ± 0.01 


< 43.0 


3 


HE 1453-0303 


0.206 


45.31 ± 0.12 


8.59 





42.54 ± 0.12 


< 43.0 


3 


SDSS 1507+4453 


0.314 


44.39 ± 0.01 


7.34 ± 0.14 


1 


42.69 ± 0.01 


< 41.66 


5 


SDSS 1507+1338 


0.322 


44.79 ± 0.01 


9.67 ± 0.02 





42.53 ± 0.01 


41.96 ± 0.08 


5 


4C 37.43 


0.371 


45.21 ± 0.01 


9.37 ± 0.11 


2 


43.18 ± 0.01 


43.02 ± 0.01 


1,2 


PG 1519+226 


0.136 


44.40 ± 0.01 


7.90 ± 0.07 





41.34 ± 0.09 


< 43.0 


1 


OR 241 


0.254 


44.86 ± 0.01 


8.44 ± 0.05 


1 


42.39 ± 0.01 


< 41.85 


2 


30R 323.1 


0.265 


45.20 ± 0.01 


9.07 ± 0.03 


2 


42.99 ± 0.01 


42.35 ± 0.13 


1, 2, 3 


40 11.50 


0.436 


44.99 ± 0.01 


8.87 ± 0.08 


2 


43.23 ± 0.01 


42.68 ± 0.03 


2 


Ton 256 


0.131 


44.51 ± 0.01 


7.99 ± 0.02 


1 


43.00 ± 0.01 


42.80 ± 0.12 


1, 3 


SDSS 1655+2146 


0.154 


44.67 ± 0.01 


8.28 ± 0.04 





43.03 ± 0.01 


43.03 ± 0.09 


3 


PG 1700+518 


0.292 


45.43 ± 0.09 


8.98 





42.63 ± 0.10 


41.81 ± 0.18 


2, 3 


30R 351 


0.372 


45.79 ± 0.01 


9.81 ± 0.04 


2 


43.50 ± 0.01 


< 42.37 


2 


n Zw 136 


0.061 


44.40 ± 0.10 


8.27 





41.95 ± 0.09 


40.90 ± 0.32 


1, 3 


PKS 2135-147 


0.200 


44.97 


9.15 


2 


43.16 


41.85 ± 0.04 


2 


HE 2152-0936 


0.192 


45.59 ± 0.10 


8.76 





42.17 + 0.11 


< 43.0 


3 


HE 2158-0107 


0.213 


44.89 ± 0.11 


8.60 





42.59 ± 0.09 


42.32 ± 0.13 


3 


HE 2158+0115 


0.160 


44.21 ± 0.10 


7.94 





42.01 ± 0.09 


< 43.0 


3 


40 31.63 


0.298 


45.99 


8.54 


1 


42.92 


< 42.08 


1, 2 


PG 2214+139 


0.067 


44.44 ± 0.11 


8.64 





41.81 ± 0.10 


< 43.0 


3 


PG 2233+134 


0.326 


45.19 ± 0.01 


8.57 ± 0.31 





42.50 ± 0.07 


< 40.58 


2 


PKS 2251+113 


0.325 


45.07 


9.15 


2 


42.98 


42.51 ± 0.04 


1, 2 


HE 2307-0254 


0.221 


44.89 ± 0.10 


8.39 





42.00 ± 0.11 


< 43.0 


3 


40 09.72 


0.433 


45.17 


9.30 


2 


42.80 


< 42.71 


2 


PKS 2349-014 


0.174 


44.77 ± 0.02 


8.81 ± 0.03 


2 


42.45 ± 0.01 


42.11 ± 0.11 


3 


HE 2353-0420 


0.229 


44.42 ± 0.10 


8.12 





42.58 ± 0.09 


42.47 ± 0.12 


3 



References. — (1) Boroson et al. 1985; (2) Stockton & MacKenty 1987; (3) Husemann et al. 2008 
(including the re-analyzed data given in Husemann et al. in prep.); (4) Husemann et al. 2010; (5) this 
work. 

'^ Radio flag. 0; RQQ, 1; FSRQ (flat-spectrum or core-dominant radio quasar), 2; SSRQ (steep-spectrum 
or lobe-dominant radio quasar). 
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included in the band, where fx.Y is the fiher transmis- 
sion at the hne Y relative to the band average and wx is 
the band width. We summarize the derived properties of 
the type-2 quasars in Table[3l All but two sources, SDSS 
1228+0050 and SDSS 1625+3106, are radio quiet (SDSS 
1625+3106 is in th e transition range of radio -quiet and 
radio-loud objects; iVillar-Martin et all 1201 ID. In addi- 
tion, we use all the type-2 quasars in iZakamska et alj 
(|2003t ) and -60,000 galaxies extracted from the SDSS 
galaxy sample at z = 0.3 - 0.6, the redshift range of the 
main sample (the type-2 quasars in Table [3]), as control 
objects. Their rest-frame color and stellar mass are cal- 
culated in the same way as for the objects in the main 
sample. 

While the above samples of type-1 and type-2 objects 
are used below to probe the different aspects of quasars, 
we should keep in mind that t he two populations can be 
intrinsically different. Elitzurl ()2012[ ) clarified this point 
by demonstrating that the two types of AGNs are drawn 
preferentially from the distribution of dust-torus covering 
factors. Since type-2 AGNs have dustier environments 
around them than type-1 counterparts, their host galax- 
ies might tend to have higher star-formation rates since 
(i) the presence of ample dust indicates active star for- 
mation producing them; and (ii) the negative feedback 
(suppression of star formation) by the AGN radiation 
might be weaker due to the dust obscuration. 

3. EMERGENCE OF EELR 

3.1. Notes on New Sample 

The Suprime-Cam [O III] line images of the newly ob- 
served quasars are presented in Figure [2] To aid visual 
inspection, the grayscale consistently represents linear 
flux levels between zero and quasar peak values while the 
white crosses indicate the PSF sizes (twice the full- widths 
at half-maximum [FWHMs]). The most spectacular fea- 
ture is found around SDSS 1507+1338. The emission 
feature, extending across ^^50 kpc along the NE - SW 
direction, is reminiscent of the ionization cones created 
by AGN radiation. SDSS 0914+0507 has apparently ex- 
tended emission too, with a knot at ^-^35 kpc away from 
the center in the NW direction. On the other hand, 
the radial profiles of SDSS 0954+0929, SDSS 1139+4020, 
and SDSS 1507+4453 are consistent with the PSFs and 
show no signs of EELR. We measure the EELR lumi- 
nosity in an annulus of inner radius 10 kpc and outer 
radius 30 kpc. The wing components of the nuclear radi- 
ations are estimated from the PSFs and subtracted. The 
incomplete transmission of the [O III] lines through the 
NAgse filter is also corrected using the filter transmission 
function. We include the uncertainties in the continuum 
subtraction and in the correction for the wing contribu- 
tion of the nuclear radiation, background sky noise, and 
the possible flux variation assumed to be 10 % due to 
the non-photometric sky condition, in the error budget 
of the derived luminosity. The results are listed in the 
last column of Table [TJ The above measurements sup- 
port our visual inspection finding detectable EELR only 
in SDSS 0914+0507 and SDSS 1507+1338. The nuclear 
[O III] luminosity (-/jn[o hi]) is also measured within a 3- 
arcsec aperture, and found to agree with the SDSS spec- 
troscopic measurements within AlogLjvfro ml = ±0.1. 

3.2. Link to AGN 
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Fig. 2. — Suprime-Cam [O III] line images of (a) SDSS 
0914-1-0507, (b) SDSS 0954-1-0929, (c) SDSS 1139+4020, (d) SDSS 
1507+4453, and (e) SDSS 1507+1338. The grayscale represents 
linear flux levels between zero and the quasar peak values, while 
the white crosses indicate twice the FWHMs of the PSFs. The 
physical scale of the panels is 50 kpc on a side. North is up, east 
to the left. 

As described in Section 1, the presence of EELR cor- 
relates with some AGN properties. Quasars with larger 
[O III] EW, broader and bumpier Balmer lines, fainter 
Fe II, and SSRQ-like radio emission have higher possi- 
bility of accompanying E ELR around them. Th e corre- 
lation is fairly strong, as IFu fc StocktonI ()2009l ) showed 
that if certain quasars are pre-selected based on their 
radio and nuclear [O III] emissions, ~50 % of them 
would turn out to have EELR as opposed to ~20 % 
with no a priori assumption. Note that the above 
emissions arise from the different parts of AGN; broad 
Balmer lines and the Fe II lines from parsec-scale BLR, 
the [O III] lines from kiloparsec-scale NLR, and ra- 
dio and EELR emissions from more extended, galaxy- 
wide scales. The correlation among them except for 
EELR has been well kn own since the seminal work by 
iBoroson fc Greenl (|1992f l. With the principal component 
analysis technique, they showed that most of the vari- 
ance in AGN spectral properties are contained in two 
sets of correlations. A strong anti-correlation between 
the measures of [O III] and Fe II is the principal com- 
ponent 1 (PC 1), or eigenvector 1, with which Balmer 
line profile is also associated. The general consensus is 
that PC 1 is driven mainly by Eddi ngton ratio (e.g., 
lBoroson|[200^ iSameshima et ai.ll2011f) . In Figure 11 we 
plot II/3 line width versus the relative Fe II strength 
fipaTT = EW(Fe II)/EW(H;3) of the SDSS quasars in 
iShen et al] (|2011f) . The both measures constitute the 
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TABLE 3 
Nuclear and host galaxy properties of type-2 quasars with EELR measurements. 







log -^N[OIII] 


Rest- frame 


Stellar mass 






Quasar 


Redshift 


(ergs-i)- 


g - i (mag)'' 


log Mi,{MQ) 


EELR'^ 


Ref. 


SDSS 0025-1040 


0.303 


42.32 


0.75 


11.0 


+1 


3 


SDSS 0123+0044 


0.399 


42.72 


1.14 


10.9 


+1 


2 


SDSS 0217-0013 


0.344 


42.34 


0.93 


11.1 


+1 


3 


SDSS 0234-0745 


0.310 


42.36 


0.91 


10.5 


-1 


3 


SDSS 0840+3838 


0.313 


42.21 


1.13 


11.4 


+1 


1 


SDSS 0849+0150 


0.376 


41.65 


1.01 


10.9 


-1 


3 


SDSS 0920+4531 


0.402 


42.63 


0.91 


11.1 


-1 


1 


SDSS 0955+0346 


0.421 


42.19 


1.24 


11.2 


+1 


3 


SDSS 1153+0326 


0.575 


43.20 


0.62 


11.2 


+1 


3 


SDSS 1228+0050 


0.575 


42.87 


0.85 


11.0 


-1 


3 


SDSS 1307-0214 


0.425 


42.51 


1.09 


11.4 


+1 


3 


SDSS 1337-0128 


0.329 


42.31 


0.90 


11.3 


+1 


3 


SDSS 1407+0217 


0.309 


42.49 


0.86 


10.5 


-1 


3 


SDSS 1413-0142 


0.380 


42.84 


1.18 


11.0 


-1 


3 


SDSS 1546-0005 


0.383 


41.77 


1.18 


11.0 


-1 


3 


SDSS 1550+3950 


0.347 


43.05 


0.60 


11.3 


+1 


1 


SDSS 1625+3106 


0.379 


42.60 


1.35 


11.1 


-1 


1 


SDSS 1726+6021 


0.333 


42.16 


0.96 


10.5 


-1 


1 


SDSS 1739+5442 


0.384 


42.01 


1.24 


11.1 


-1 


1 


SDSS 2358-0009 


0.402 


42.91 


0.88 


11.1 


+1 


3 



References. — (1) Humphrey et al. 2010; (2) Villar-Martin et al. 2010; (3) Villar-Martm 
et al. 2011. 
'^ The measurement error is <0.01. 

Typical error of a few tens of magnitude is expected, including the uncertainties in the 
photometry and fc-correction. 
'^ EELR detection flag (+1: detected, —1: non-detected). 



(Fe II is emitted from the outermost part of 
, IMatsuoka et al.l 120081 ) as the viewing angle 



four-dimensional eigenvector 1 (4DE1; other measures 
are the soft X-ray phot on index and C I V A15 49 line pro- 
file) space proposed bv lSulentic et al.l ()2000t ). The anti- 
correlation between them is evident in Figure [31 which 
forms a sequence of Eddington ratio. The data scatter re- 
flects the orientation effects as well as the intrinsic variety 
of quasars, since one would observe larger line FWHM 
and i?Feii 
BLR; e.g. 

relative to the dust-torus plane decreases. On the other 
hand, PC 2 links optical luminosity and He II A4686 
line emission. Its physical origin is thought to be mass- 
accretion rate to SMBHs because of its strong correla- 
tion with optical luminosity. Therefore one would expect 
that the ratio of PC 1 and PC 2 is proportional to SMBH 
mass, whi ch was actua lly proved by Boros on (2002). Fur- 
thermore, iLaoii (|2000D found that the radio loudness is 
strongly related to SMBH mass as radio-loud/radio-quiet 
quasars are associated with most/least massive SMBHs. 
The above arguments suggest that most of the AGN 
properties can be understood on the PC 1 - PC 2 plane. 
We plot Eddington ratio and continuum luminosity of 
the compiled t ype-l sampl e , as w ell as the whole SDSS 
quasars in the iShen et al] ()2011l) catalog, in Figure [H 
First of all, we remark that the objects with EELR 
measurements cover the distribution of the whole SDSS 
quasars fairly uniformly on this plane. It was not an- 
ticipated since the existing EELR observations are sub- 
ject to strong selection biases individually. We see the 
clear separatioi i of radio-qu iet and radio-loud sources as 
pointed out bv iLaorl (|2000f ). the latter being associated 
with most massive SMBHs. In addition, SSRQs are ob- 
served to have higher values of Eddington ratio or SMBH 
mass than FSRQs. We also find some correlations con- 
taining EELR. They are shown more clearly in Figure 
[5] where the fraction of EELR detection is plotted as a 
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Fig. 3.— H/3 line width (FWHM) versus the relative Fe II 
strength Rpell = EW(Fe II)/EW(H/3) of the SDSS quasars. The 
red, green, and black dots represent the objects with high (—0.8 < 
log [ibolZ-f'Edd]), middle (-1.2 < log [Lbol/iEdd] < -0.8), and 
low (log [ibol/^Edd] < —1.8) values of Eddington ratio, respec- 
tively. The typical error of each data point is indicated at the 
top-right corner. 

function of the nuclear properties. The EELR fraction 
is strongly anti-correlated with Eddington ratio, and its 
weak correlation with SMBH mass may also be present. 
The former relation is most likely the underlying basis 
of the correlation between EELR and PC 1. While no 
clear trend is found with respect to radio emission, we 
do find that EELR luminosity is systematically higher in 
SSRQs than in the rest of the sample; the mean luminos- 
ity is Le[oiii] = 10*^* and lO^^-^ erg s^^ for the SSRQs 
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and RQQs, respectively. The difference would be even 
larger when we consider the fact t hat the EELR luminos- 
ity of most SSRQs measured by iStockton fc MacKentvl 
(|1987| ) is underestimated since the emission in a cen- 
tral 10-kpc aperture is excluded. It explains why the 
EELR fraction around SSRQs is much higher than those 
around the other popu l ations in the pioneering work of 
IStockton fc MacKentvl (|1987D . 

One of the highest dependency of EELR is found on 
nuclear [O III] luminosity, ^^[0111]- In Figure [6] we 
show Eddington ratio versus -^n[oiii] of the same sam- 
ples. It is clearly seen that the quasars with EELR 
measurements are biased to high Ln[oiii] compared to 
the whole SDSS quasars. Here we find that 1ow-Ln[oiii] 
quasars tend to lack EELR: only 3 out of 20 objects with 
-^N[oiii] < lO"'^'^ erg s~^ have detectable, faint EELR, 
while the EELR fraction is > 50 % at the higher iN[oiii]- 
This trend was reported previous ly for smaller samples 
(e.g.. IStockton fc MacKentvlll987[) . Figure H] also tells us 
that iN[oiii] does not strongly correlate with Eddington 
ratio , although it correlates with PC 1 (jBoroson fc GreenI 
119921) . No clear correlation between the two measures 
indicates that PC 1 is not solely driven by Eddington 
ratio - there should be another physical mechanism(s) 
that links nuclear [O III] luminosity with other PC 1 
constituents. 

What can we learn from these correlations? While the 
clear anti-correlation between EELR and Eddington ra- 
tio is observed, it does not seem to reflect their direct re- 
lationship through the AGN radiation since (i) the EELR 
fraction is highest for the AGNs with most inefficient gas 
accretion to SMBHs; (ii) nuclear [O III] luminosity pro- 
duced similarly by the AGN radiation shows no clear cor- 
relation with Eddington ratio; (iii) a positive correlation 
between the EELR fraction and continuum luminosity is 
not observed. Therefore it is more reasonable to consider 
the presence of an independent mechanism that leads to 
the high/low EELR fraction and low/high Eddington ra- 
tio at the same time. The dependence of EELR on radio 
emission and (possibly) SMBH mass may also be related 
to this hidden mechanism. We will come back to this 
issue further below. The lack of simple relation between 
EELR and conti nuum luminosity (acc retion rate) was 
also reported bv iFu fc Stockton! ()2009( ). who suggested 
that the correlation between nuclear and extended [O III] 
emissions, which is also evident in our larger sample, may 
be the sequence of the available gas amount rather than 
AGN emission characteristics. The hypothesis that the 
gas availability is the primary determinant of EELR cre- 
ation is supported by the probe of host galaxies in the 
following section. 

3.3. Link to Host Galaxies 

If EELR represents a certain stage of galaxy evolution, 
as the AGN feedback scenario implies, their presence may 
be related to stellar-population properties of host galax- 
ies. Here an investigation of such a relation is presented 
for the first time, thanks to the recent measurements 
of EELR around SPSS type-2 quasars at z = 0.3 - 0. 6 
(IHumohrev et aLllMl rVinar-Martin"eraIll20Toll20TTI) . 
In Figure [71 we plot stellar mass versus rest-frame g — i 
color of the type-2 quasars as well as the control samples 
defined in Section \T2\ Most SDSS galaxies at these rel- 



atively high redshifts are luminous red galaxies (LRGs), 
forming a clump at the massive red end on this color 
- stellar mass diagram (we could see the red sequence 
rather than a clump if the SDSS observations were deep 
enough to detect fainter galaxies ). The type-2 quasars 
are bluer than the LRGs (see also. lZakamska et al.ll2003l ) 
but are redder than star-forming galaxies which form the 
"blue cloud" around g — i ^ 0.6 in the local universe (e.g., 
Gavazzi et al.ll2010fl. It is consistent with the result of 



Schawinski et al.l (|2010( ) who found that local AGNs re- 



side in the "green valley" between the blue cloud and 
the red sequence on the color - stellar mass diagram. 
Interestingly, we find that the EELR quasars are asso- 
ciated with bluer and more massive galaxies compared 
to the non-EELR counterparts. In order to show this 
trend more clearly, we define a 'massive-blue'-ness indi- 
cator, Imb = log Mj, — [g — i), which increases diagonally 
toward the bottom-right corner of Figure [7] as shown by 
the arrow. The number distribution of the type-2 quasars 
with EELR measurements as a function of Imb is plot- 
ted in Figure |8l It shows a clear, monotonic increase in 
the fraction of EELR quasars as X^b increases: only 1/8 
of the lowest-Imb quasars has EELR, while the fraction 
rises to 7/8 in the highcst-I,nb counterparts. 

What is the origin of this dichotomy? Blue colors point 
to current or recent star-formation activity, which in turn 
indicates the presence of a certain amount of gas in the 
galaxies. Based on the velocity structure and mass of 
EELR gas around mostly SSRQs, as well as the correla- 
tion b etween EELR and radio emission, iFu fc Stockton! 
(|2009l ) suggested that EELR is formed in galaxy minor 
merger which is accompanied by production of a radio jet 
associated to AGN. A s imilar conclu sion was reached by 
'Huse mann et al.l (|2010f ) for a RQQ. iVillar-Martin et al.l 
(2011') found that the majority of the type-2, radio-quiet 
EELR quasars show the signs of galaxy merger or inter- 
action while none of the non-EELR counterparts show 
such signs. These observations suggest that the gas 
may be brought in by the merging companions. The 
merger/interaction possibly activates the observed star 
formation as well as AGN that ionizes the gas and creates 
EELR, since the interstellar gas of merging galaxies can 
be effectively brought into the nucle ar regions with rapid 
loss o f the angular momentum (e.g.. lBarnes fc Hernguisti 
119960 . The fact that EELR is preferentially found in 
gas-rich, massive blue galaxies is consistent with both the 
above scenario and the previous arguments that the pres- 
ence of EELR is regulated by the amount of available gas. 
Furthermore, massive galaxies are kno wn to be the dom- 
inant hosts of r adio-loud AGNs (e.g.. iFlovd et al.l!200l 
iBest et al]|2005f ). It naturally explains the EELR - radio 
emission connection described above, and can also pro- 
vide the " hidden" mechanism behind the anti-correlation 
between EELR and Eddington ratio since the majority of 
radio-loud AGNs are powered by inefficiently-accreting, 
most massive SMBHs (see Figure HI). 

The recent kinematic observations have revealed the 
presence of g alaxy- wide, disturbed gas in quasar host 
galaxies (e.g.. lFu fc Stockton 2009; Gr eene et a l. 2011). 
While only a small fraction of the gas is generally found 
to have high enou gh velocity to escape the galaxies, 
iGreene et al.l (|201in showed that these estimates are in 
fact lower limits and may be consistent with expectations 
of recent AGN feedback models. If the feedback scenario 



Matsuoka 



44.0 - 



o 
o 

f 44.5 



en 
O 



g 45.0 

c 



^ 45.5 

c 

c 

o 

O 



46.0 - 



PC2 
A 








MgH.^ const 



® 


EELR 


• 


no EELR 




SDSS 


• 


RQQ 


• 


FSRQ 


• 


SSRQ 



J I I I • I -i I L 



-> PCI 

_I I I I I I I I I I L 



0.5 



0.0 -0.5 -1.0 -1.5 -2.0 

Eddington ratio (log L^jL^^^) 



Fig. 4. — Eddington ratio versus continuum luminosity of the compiled type-1 sample (filled circles) and tfie whole SDSS quasars (dots). 
RQQs, FSRQs, and SSRQs in the compiled sample are represented by black, green, and red colors, respectively. The quasars with detectable 



EELR are marked with the larger circles whose sizes are proportional to 
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PC 1 and PC 2 increase along the horizontal and 



vertical axes on this plane, while SMBH mass increases from the top-left to the bottom-right corner, 
log Mbh = 8-0> 8.5, 9.0, 9.5, are shown by the dotted lines. 



The lines of constant SMBH mass. 



is correct, the gas responsible for EELR and star forma- 
tion (hence blue colors of galaxies) would soon be swept 
away by the AGN-driven, galactic winds. In this regard, 
it may be noteworthy that none of the most luminous (log 
AiAsioo ^ 45.5) type-1 quasars in our sample has EELR 
as shown in Figure SI It is tempting to raise the possi- 
bility of the gas stripping, but we have to be aware that 
(i) the relatively-high Eddington ratio for these objects 
may be responsible for the deficit of EELR; (ii) the bright 
central radiation can prevent us from detecting EELR in 
the most luminous sources. More data is needed for fur- 
ther inspection on this issue. While the process discussed 
above corresponds to the so-called " quasar- mode" feed- 
back, it is interesting to point out that a clear interaction 
between radio jets and EELR clo uds with high velocit y 
dispersion is olDserved in SSRQs ()Fu fc Stocktonll2009l ). 



which invokes the another , "radio mode" of the feedback. 
In the nearby universe, iKeel et all ()2012D investigated 
the properties of Seyfert galaxies with EELR, dominated 
by type-2 objects. While their morphology supports the 
merger origin of EELR, they have relatively quiescent 
kinematics without a sign of high- velocity gas leaving the 
galaxies seen in higher-luminosity AGNs. It may sim- 
ply indicate that creation of the powerful galactic winds 
requires quasar-class energy input into the interstellar 
medium, hence the relation of EELR to AGNs may be 
different at high and low luminosity. 

4. SUMMARY AND CONCLUSIONS 

A comprehensive analysis of EELR phenomenon 
around quasars is presented in this paper. We 
compile the past EELR measurements for type- 
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Fig. 5. — Fraction of EELR detection as a function of Eddlngton 
ratio (panel a), SMBH mass (b), continuum luminosity (c), and 
radio-emission property (d). The fractions are calculated in the 
bins represented by the horizontal bars. The number of objects In 
each bin is also indicated. 



1 (Boroson et al.' 1981 iStockton fc Ma cKentv '1987"; 
Husemann ct al. 2008) and t ype- 2 ^Humphrey ct al. 
20101 : iVillar-Martin et al.ll2010l feoill ) quasars, and com- 
bine them with the new obseryation of fiye quasars 
at z '^ 0.3. The new obseryation was carried 
out using the Subaru/Suprime-Cam with the narrow- 
band NAgse filter (corresponding to the redshifted 
[O III] A5007), and led to the discoyery of EELR 
around two sources, SDSS 091401.75-^050750.6 and 
SDSS 150752.66-^133844.5. The observation effectiyely 
complements the existing measurements on the AGN 
principal-component (PC) plane. The properties of asso- 
ciated AGN (SMBH mass, continuum and line luminos- 
ity, and radio characteristics) and host galaxies (stellar 
mass and rest-frame color) are collected and calculated 
from the SDSS and other obseryations. 

We find that EELR anti-correlates with Eddlngton 
ratio, which is most likely the underlying basis of the 
known correlations between EELR and the PC 1 (eigen- 
ycctor 1) constituents. EELR and nuclear [O III] emis- 
sions are also highly correlated. We argue that the pri- 
mary determinant of EELR creation is the gas availabil- 
ity, rather than AGN em ission character istics. The re- 
cent o bseryations (iFu fc Stockton 2009; H usemann et al.l 
I2010t IVillar-Martin et al.l 120111 ) revealed that a signif- 
icant fraction of EELR galaxies show the signs of re- 
cent minor merger, which may have brought in the gas 
and activated star formation as well as AGN that creates 
EELR. In line with the above arguments, we find that 
EELR is preferentially associated with gas-rich, massive 
blue galaxies. It may explain the whole set of observed 
correlations containing EELR, since massive galaxies are 



the dominant hosts of radio-loud AGNs characterized 
with low Eddlngton ratio and high SMBH mass. 

The hierarchical structure formation models based on 
the ACDM theory predict that galaxies are assembled 
through a sequence of mergers of smaller building blocks. 
In Figure [9] we show a schematic view of the color - stel- 
lar mass diagram, or equivalently the color - magnitude 
diagram (CMD), of observed galaxies. Once a galaxy 
is grown up to the massive tip of the blue cloud (rep- 
resented by the seven-rayed star in the figure), it has 
two ways to evolve further. If no merger happens sub- 
sequently for sufficient time, the galaxy would migrate 
into the low-mass end of the red sequence as the star for- 
mation fades out. It can then move up the red sequence 
toward the massive en d through dry me rger with gas- 
poor galaxies (see also. lFaber et al.ll2007[ ). On the other 
hand, wet merger of the galaxy with gas-rich galaxies at 
any stage of the above evolution could push the galaxy 
into the green valley where AGN fraction is the high- 
est. If the galaxy (and hence the central SMBH, assum- 
ing the positive bulge mass - SMBH mass correlation) 
is massive enough and the merging events involve am- 
ple supply of gas for triggering vigorous star formation, 
the galaxy would be driven to the massive blue realm 
at the extreme corner of the valley. This event would 
also activate AGN, which ionizes the gas and creates lu- 
minous EELR. The galaxy-wide energy injection by the 
AGN radiation might then blow away the gas and quench 
the star formation, pushing the galaxy back to the red 
sequence. While this last argument is less certain due 
to the lack of evidence about the energy injection pro- 
cess, recent observations have started to reveal the pres- 
ence of AGN-driven, galaxy-scale outflows which may be 
quenching the star formation there, fe.g.. 'Greene et ed\ 
2011; Rupke & Veilleux 2011; Cano-Diaz ct al. 20U^! 

The above AGN feedback scenario gives a solution to 
the major problems of the classical hierarchical galaxy 
formation models, i.e., the overproduction of massive 
galaxies and the inverted color - magnitude - morphology 
relation (see Section 1) by suppressing star formation in 
massive galaxies. However, many details of the process 
are still unknown and to be debated. What is observa- 
tionally clear is that EELR is found at the exact region on 
the CMD where the AGN feedback is expected to work 
- the massive blue corner that seems to be the forbidden 
area for observed galaxies. It provides a piece of evi- 
dence for the presence of the AGN feedback, which may 
be playing a leading role in the co-evolution of galaxies 
and central SMBHs. 
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